Correlated magnetic states in domain and grain boundaries in graphene 
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A b initio calculations indicate that while the electronic states introduced by grain boundaries 
in graphene are only partially confined to the defect core, a domain boundary introduces states 
near the Fermi level that are very strongly confined to the core of the defect, and that display a 
ferromagnetic ground state. The domain boundary is fully immersed within the graphene matrix, 
hence this magnetic state is protected from reconstruction effects that have hampered experimental 
detection in the case of ribbon edge states. Furthermore, our calculations suggest that charge 
transfer between one-dimensional extended defects and the bulk in graphene is short ranged for 
both grain and domain boundaries. 

PACS numbers: 73.22.-f, 73.20.Hb, 71.55.-i 



Controlling electronic transport and tailoring magnetic 
states at the nanoscale in graphene rank among the main 
issues related to the prospective application of this ma- 
terial in nanoelectronics and spintronics. In particular, 
electronic and magnetic states of extended line defects 
in graphene have been considered as possible conducting 
one-dimensional (ID) electronic channels and platforms 
for tailored spin states for spintronic applications, p], Q 
Two recent developments highlight the focus on extended 
line defects: (i) the recognition that mass-scale produc- 
tion of graphene should inevitably lead to a polycrys- 
talline material, containing one-dimensional grain bound- 
aries (GB) |3|; and (ii) the recent theoretical predic- 
tion [4] and experimental realization of domain bound- 
aries (DB) in graphene by controlled deposition of the 




FIG. 1: Geometries of extended one-dimensional defects in 
graphene. Carbon (oxygen) atoms are shown as grey (black) 
circles, (a) DB558: a domain boundary, (b) DB558+0: oxi- 
dized form with one O atom per defect unit, (c) DB558+20: 
oxidized form with two O atoms per defect unit, (d) GB: a 
large- angle tilt grain boundary, (e) GB+O: oxidized form of 
the GB. (f) Schematic view of the first quadrant of the Bril- 
louin zone associated to the defect supercells. High-symmetry 
points and Cartesian axes are indicated. 



material on metallic substrates [2[. In the case of GB's, 
a large volume of works has accumulated in the last few 
years [l-[s| , while for the DB produced in Lahiri et al. ex- 
periment, valley- filter properties [9] and its effects on the 
magnetic edge states of a graphene ribbon [To| have been 
theoretically investigated. Further, when chemical reduc- 
tion of a readily available material, such as graphene ox- 
ide, is used as a viable route for mass-scale graphene syn- 
thesis [UE2, the presence of residual functional groups 
should affect the material electronic and magnetic prop- 
erties. 

In this context, the nature of the electronic states in- 
troduced by such extended ID defects in graphene is a 
topic that deserves close inspection. More specifically, 
whether GB's and DB's act as quasi- ID conducting chan- 
nels immersed in the bulk of graphene is the question we 
seek to address in this work. We also consider the is- 
sue of self-doping in graphene induced by the presence 
of such extended ID defects [Hj], that occurs when the 
line defect attracts charge carriers, resulting in charged 
defective lines surrounded by a doped graphene matrix. 

Our calculations indicate that while the electronic 
states introduced by a model GB structure in graphene 
hybridize with the bulk states and are only partially con- 
fined to the defect core, the DB defect introduces a sharp 
resonance in the density of states (DOS) of graphene, 
that lies just above the Fermi level in the neutral sys- 
tem, and is associated to electronic states that are very 
strongly confined to the core of the defect. Our re- 
sults suggest that, when a graphene sample containing 
a DB is doped and these quasi- ID states are populated, 
a ferromagnetic state is realized, which is confined to 
zigzag chains along the defect core that are fully im- 
mersed within a bulk graphene matrix. The quantum 
confinement induced by the presence of a domain bound- 
ary leads to an enhancement of the Coulomb interactions 
and stronger electron-electron correlations. Because of 
their ID nature, these correlated states do not show long 
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range order. Instead, they present power law or algebraic 
correlation functions. Furthermore, given that these ID 
states are immersed in a metallic environment, they are 
unique examples of open Luttinger liquids [13 ]. 

Moreover, we find that charge transfer between bulk 
graphene and the line defects is essentially local, with 
charge redistribution taking place within a region of 
~3-5 A from the geometric center of the ID defects. 
Both this charge-transfer region and the surrounding 
bulk regions remain essentially neutral. The absence of 
a charge monopole moment means that these ID defects 
should act as weak charge carrier scatterers, in agreement 
with recent transport experiments in graphene grown by 
chemical vapor deposition (CVD) [15j. 

In this study, we examine the nature of the electronic 
states introduced by the two aforementioned types of 
extended ID defects in graphene: (i) large-angle tilt 
grain boundaries p, Q , expected to occur commonly in 
mass-scale production of graphene; (ii) domain bound- 
aries produced by deposition of a graphene layer on a 
Nickel substrate [2|, |4j. For both ID defects, we also con- 
sider oxidized forms, with oxygen atoms bound to the 
atoms at the core of the defect. Figure QJa) shows the 
atomistic model for the domain boundary defect, which 
consists of periodic units composed of one octagon and 
two side-sharing pentagons along the defect line, which 
we label DB558 in the following discussion. The oxi- 
dized forms with one (DB558+0) and two oxygen atoms 
(DB558+20) per defect unit are shown respectively in 
Figs.QJb) and (c). Figure QJd) shows the GB model pro- 
posed in Ref. 6, while Fig. [TJe) shows the fully unzipped 
oxidized form of this defect (GB+O). 

Our calculations are performed in the framework of 
Kohn-Sham density functional theory (DFT), within 
the generalized-gradient approximation (GGA) [l6| and 
norm-conserving pseudopotentials in the Kleinman- 
Bylander factorized form [17j. We use the LCAO method 
implemented in the SIESTA code [Hj , with a double-zeta 
basis set plus polarization orbitals. A Mulliken charge 
partition is employed for the analysis of the charge- 
transfer between bulk and defects. Supercells of 66 (84) 
atoms are employed for the DB558 (GB) calculations. 
We performed convergence tests to ensure that, at these 
cell sizes, our results are not affected by interaction be- 
tween the defects and their periodic images. 

In the case of GB's in graphene, the electronic struc- 
ture has been amply discussed in Ref. 5, where the oc- 
currence of an anisotropic Dirac cone in a Brillouin-zone 
(BZ) line along the GB direction was analyzed. In Ref.y, 
the GB electronic states near the Fermi level (FL) were 
found to disperse in all directions, indicating hybridiza- 
tion with bulk states. Further indication of such hy- 
bridized nature of GB states is provided by the analysis 
of the full band structure and density of states (DOS) 
of the GB supercell. Figure [2ja) shows the band struc- 
ture along the BZ lines shown in Fig. DJf). Note that 
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FIG. 2: (Color online) (a) Band- structure for the GB along 
the Brillouin-zone lines in Fig. 1(f). (b) DOS for the GB. (c) 
Band-structure for the DB558. (d) DOS for the DB558. 

the states near the Dirac point show sizeable dispersions 
in directions perpendicular to the defect (Y-M and T-X 
lines). The contributions of the orbitals centered on the 
atoms at the GB core to the resonances introduced in 
the DOS of graphene are shown in Figure EJb). For the 
purpose of our analysis, we define the core of both the 
GB and DB558 defects as composed by the ten atoms 
that form the topological defects along the defect line, as 
indicated in Fig. [TJ In the DOS plots, we refer to atoms 
1 and 2, that are placed along the line at the geometric 
centers of both defects, as core center, while the remain- 
ing eight core atoms are referred to as core perimeter. 
In the case of the GB, the contribution of the ten core 
atoms to the four DOS peaks near the FL ranges from 
41% to 57%. These results will be contrasted with the 
case of the DB558 in the following. 

The band structure and the DOS for the DB558 are 
shown in Fig. [2fc) and (d), respectively. In Fig. [2fc) 
we also include the band structure of a 64-atom bulk 
supercell, obtained by removing core atoms 1 and 2 from 
the geometry in Fig. QJa). Given the periodicity of both 
defect and bulk cells along the DB558 direction (a zigzag 
direction of the bulk matrix), the Dirac point of the bulk 
cell folds onto the point at 2tt/3£, where i is the period of 
the DB558 [19[. Note that, to a large degree, the changes 
in the DB558 band structure with respect to the bulk 
one are concentrated in a region of ±0.5eF from the FL, 
and that within ~0.2 eV above (below) the FL we find 
bands with flat sections along the T— Y (L-X) direction 
(both parallel to the defect line). These states show little 
or no dispersion along the T-X and Y-M lines that are 
perpendicular to the defect line. The flat-band character 
of the states above the FL, near the zone center, leads 
to a ferromagnetic state for a bulk ID defect in graphene 
consisting entirely of carbon atoms, as discussed below. 

Considering now the DOS for the DB558 in Fig.^d), 
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FIG. 3: (Color online) Isosurface of spin polarization density 
for domain-boundary magnetic state. 



we find two sharp resonances within ~ ±0.1 eV from the 
FL, associated to the ID defect. These resonances reflect 
the presence of extended van Hove regions in the in the 
bandstructure, as seen in Fig. 2(c). We observe that 
in one feature the DOS of the DB558 differs markedly 
from the GB one: the peak just above the Fermi level 
shows a very strong concentration on the core atoms, 
with ^80% of the total DOS concentrated on the orbitals 
centered on the perimeter core atoms numbered 3, 4, 5, 
and 6, along the zigzag chains that bond to the dimers 
(atoms 1 and 2) at the core center in Fig. QJb). Adding 
the contributions of the atoms from the same sublattice 
as atoms 3-6, in the zigzag chains nearest to the core 
on both sides of the defect line, we already account for 
~96% of the total DOS for this peak. A strongly one- 
dimensionally confined empty state along the defect core 
is thus a characteristic of the DB558. In our calculations, 
however, these states are empty in the neutral system, 
hence no spin polarization is induced. 

In order to investigate the formation of magnetic states 
in gate-doped versions of this system, we added one elec- 
tron per supercell, a doping concentration that raises the 
FL by about ~0.07 eV. A spin-polarized calculation sta- 
bilizes a ferromagnetic state, with a magnetic moment of 
0.52 \±b per defect unit, and a formation energy that is 
~40 meV per defect unit lower in energy than the unpo- 
larized state. Figure [3] shows a representative isosurface 
of the difference between majority and minority spin den- 
sities. The spin density is concentrated on the sublattice 
of the atoms 3-6, with a negligible contribution from the 
dimer atoms themselves and also from the atoms near 
the core that belong to the other sublattice. We see here 
a manifestation of a magnetic state along a line defect 
that is fully immersed within the bulk of graphene. The 
fact that these states have a negligible contribution along 
the dimer atoms is perhaps the origin of the topological 
disruption in the electronic states that leads to the local- 
ization of the related 7r-orbital bands and hence to the 
magnetic state. We speculate that such DB magnetic 
states should be more easily detectable experimentally 




FIG. 4: (Color online) (a) Band-structure for the GB+O. (b) 
Band-structure for the DB558+20. (c) Band-structure for the 
DB558+0. (d) Band-structure for the DB558+0, including 
spin polarization. Black (red) symbols show majority (minor- 
ity) spin bands. 



than those predicted to occur along the edges in graphene 
ribbons [19l-l22| , since the zigzag chains along the DB are 
protected from reconstruction. 

To address the effect of chemical doping on the elec- 
tronic states of the ID defects, we considered the oxi- 
dized forms shown in Fig. HJ The corresponding band 
structures are shown in Fig. HI For both the GB-j-0 and 
DB55+20, where oxygen-induced unzipping of the C-C 
bonds takes place, we observe the opening of small gaps 
of 55 meV and 190 meV, respectively. In neither of these 
two cases we find magnetic ground states in the neutral 
system. An interesting case is the DB558+0 system, 
with one oxygen atom per defect cell, in which the oxy- 
gen atom bonds with the two carbon atoms at the defect 
core in the bridge position. In this case, as shown in 
Fig. E^c), the occupied band along the T— Y line that 
is ~0.1 eV below the FL in the pure DB558 system is 
now deeper in energy. More importantly, a significant 
change is observed in the flat band of confined ID states 
near T that now crosses the FL. As a result, a magnetic 
ground state is stabilized, with a net magnetic moment of 
l.lfiB per unit cell, as shown in Fig. H^d), where major- 
ity and minority spin bands are plotted. Note that spin- 
polarization is strongly localized on the ID defect, with 
significant exchange splitting restricted to the quasi- ID 
bands crossing the FL. This ferromagnetic ground- state 
is 43 meV (per defect unit) lower in energy than the 
non-polarized state. This value is very close to that for 
the pure DB558 system with an extra electron per cell. 
Hence, chemical doping offers another way of stabilizing 
this ferromagnetic state. The ferromagnetic instability 
of these flat-band states is further confirmed by a calcu- 
lation we performed for the DB558+20 system with an 
extra charge added to the defect supercell. Again, the FL 
shifts into the band of quasi- ID states, and the system 
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FIG. 5: (Color online) (a) Charge density distribution for the 
DB558 (blue circles) and the GB (red squares), (b) Integrated 
charge J(r) (see text) for the DB558 and the GB. 

stabilizes a ferromagnetic state. Just as in the case with- 
out O doping, the presence of extended van Hove regions 
is also clearly seen in the band structures in Fig. 2J 

We turn our attention now to the charge transfer be- 
tween the ID defects and the surrounding graphene bulk, 
since this may give important qualitative information 
about the nature of the scattering of charge carriers by 
such ID defects. In Fig. [5j we show the profile of the 
charge distribution around the ID defects, and the inte- 
grated charge per atom, both as a function of the distance 
to the geometric center of the defect, for the DB558 and 
the GB. The charge distribution is shown as linear charge 
densities representing the net charge summed over atoms 
at the same distance from the ID defect, divided by the 
period of the defect unit. The integrated charge is the 
integral of this linear charge density profile, computed as 
total net charge per atom, and defined as follows: 

R 

HR)= N(Bj E A ^)x^; (1) 

where A(r^) is the linear charge density for the line of 
atoms at a distance n from the defect, £d is the period of 
the ID defect, and N (R) is the number of atoms summed 
over all lines with T{ < R. 

Figure 5 shows that the charge redistribution around 
the ID defects is non-monotonic, with the linear charge 
distributions alternating in sign on both sides of the de- 
fect, as seen in Fig. [5fa). Note also that charge trans- 
fer between the ID defects and the graphene matrix oc- 
curs in a range of ~3-5 A from the defect line, with the 
graphene matrix becoming neutral beyond this range, 
as shown in the plots for I(R) in Fig. [5jb). Since the 
charge distribution region recovers neutrality within ~3- 
5 A from the defect center, we expect no long range 
Coulomb scattering of charge carriers from both ID de- 
fects. In our calculations, oxidation of the ID defects is 
found not to affect this localized character of the charge 
balance between the defects and the graphene bulk. Such 
extended ID defects should then play a minor role as a 
source of carrier scattering in graphene. 

In summary, we find that the electronic states intro- 
duced by grain boundaries in graphene are only partially 



confined to the defect core, while a domain boundary in- 
troduces unoccupied electronic states near the Fermi level 
that are very strongly confined to the core of the defect, 
and that, when populated by doping, display a ferromag- 
netic ground state in a ID defect that is fully contained 
within the bulk matrix and that consists entirely of car- 
bon atoms. Being fully bulk-immersed, this ferromag- 
netic state is protected from reconstruction and should 
be more easily detectable experimentally than those pre- 
dicted to exist along the edges of graphene ribbons. Fur- 
thermore, our calculations indicate that charge transfer 
between bulk graphene and both types of extended ID 
defects is strongly localized, with charge redistribution 
confined to ~3-5 A from the geometric center of the ID 
defect, implying that these ID defects should act as weak 
charge scatterers in graphene. 
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